Influence of gauge artifact on adiabatic and entropy perturbations during inflation 
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In recent publications we proposed one way to calculate gauge-invariant variables in the local observable 
universe, which is limited to a portion of the whole universe. To provide a theoretical prediction of the observable 
fluctuations, we need to preserve the gauge-invariance in the local universe, which is referred to the genuine 
gauge-invariance. The importance of the genuine gauge-invariance is highlighted in the study of primordial 
fluctuations in the infrared (IR) limit. The stability against IR loop corrections to primordial fluctuations is 
guaranteed in requesting the genuine gauge-invariance. The genuine gauge-invariance also gives the impact on 
the detectable fluctuations. We showed that at observable scales the bi-spectrum calculated in the conventional 
perturbation theory vanishes in the squeezed limit, if we request the genuine gauge invariance. This indicates 
that the conventional bi-spectrum is, in this limit, dominated by a gauge artifact, which cannot be observed. 
These studies have been elaborated in single-field models of inflation. In this paper we generalize our argument 
to multi-field models of inflation, where, in addition to the adiabatic field, the entropy field can participate in the 
generation of primordial fluctuations. We will find that the entropy field can generate the observable fluctuations, 
which cannot be eliminated by gauge transformations in the local universe. 



I. INTRODUCTION 

The conventional cosmological perturbation theory has 
been established based on the assumption that we know the 
whole spatial region of the universe with infinite volume. We 
should, however, recognize that this does not meet the case 
in the actual observations, because the observable portion of 
the universe is limited. In our recent works (Hl-01, we pointed 
out the necessity of distinguishing the gauge invariance in the 
whole universe with infinite volume from that in the local uni- 
verse with finite volume. To preserve the gauge-invariance in 
the whole universe, it is sufficient to request the invariance un- 
der normalizable gauge transformations, which become regu- 
lar at spatial infinity. However, in the local universe, where 
we need not to concern about the regularity at infinity, it is 
necessary to request the invariance under both normalizable 
and non-normalizable gauge transformations. We discrimi- 
nate gauge-invariant quantities that can be constructed in our 
local observable universe, referring them as genuine gauge- 
invariant variables. The observable fluctuations should be 
such a genuinely gauge-invariant variable. 

The genuine gauge-invariant perturbation has a large im- 
pact namely on the infrared (IR) modes of fluctuations. The 
adiabatic vacuum, which yields the scale-invariant spectrum, 
is supposed to be a natural vacuum in the inflationary uni- 
verse. However, once the interaction turns on, it is not mani- 
fest whether the adiabatic vacuum is stable or unstable against 
the IR contributions in loop corrections JH-H^]. (See also 
the recent discussions in Refs. ll27l - [35ll .) In our previous 
works [2, 3], we showed that in single field models of inflation 
the IR divergence is an unphysical artifact, which disappears 
in requesting the genuine gauge-invariance. Namely, initial 
quantum states are requested to satisfy several conditions in 
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order to respect the gauge invariance in the local universe. 
The importance of the genuine gauge invariance is recognized 
as well in the study of the primordial non-Gaussianity fl. 
It is remarkable that the tree-level bi-spectrum calculated in 
the conventional perturbation theory vanishes in the squeezed 
limit under the request of the genuine gauge invariance. This 
fact emphasizes the importance to investigate primordial fluc- 
tuations based on the genuine gauge invariant perturbation 
theory to yield the theoretical prediction of fluctuations, which 
are to be compared with the observations. 

These arguments have been so far elaborated in single field 
models of inflation. In this paper, we extend our argument to 
multi-field models. In multi-field models, the issue of gauge- 
invariance becomes more delicate because of the presence of 
the entropy field. If the entropy field is massless, the loop cor- 
rection of the entropy field diverges due to IR contributions. 
In contrast to the IR divergence from the adiabatic field, the IR 
divergence from the entropy field is conceived to be irrelevant 
to the presence of the gauge degrees of freedom J36|. The 
purpose of this paper is not to provide ways of regularization, 
but to provide one way to realize the genuine gauge invari- 
ance in the presence of the entropy field. For this purpose, 
we need to discriminate between the IR divergences from the 
gauge effects and those from different origins. In this paper 
we consider two-field models of inflation, but an extension 
to inflation models with more than two fields would proceed 
straightforwardly. 

Our paper is organized as follows. In Sec. HU we give the 
setup of our problem and after that we briefly review a method 
to calculate genuinely gauge-invariant quantities. In Sec. [HI] 
we study the implications of the genuine gauge-invariant per- 
turbation in the two-field models with the pure adiabatic and 
entropy fields. After we derive the gauge-invariance condi- 
tion, we calculate the genuinely gauge-invariant bi-spectrum 
in this model. In Sec. IIV1 we extend our arguments to more 
general two-field models where the adiabatic and entropy 
fields are coupled even at liner order. Our results are sum- 
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marized in Sec. [VI 



B. Gauge-invariant operator 



II. GENUINE GAUGE-INVARIANT PERTURBATIONS 

In this section, we first give the set up of the problem. Af- 
ter that, we provide one way to perform the genuinely gauge- 
invariant perturbation theory, following the discussion in our 
previous works |0, . 



A. Basic equations 

We consider the two-field inflation models with the stan- 
dard kinetic term whose action take the form 

M 2 r 

S = ^T] V=g[R- gT^i,* ~ 2^(01, <h)] d 4 x, 

(2.1) 

where M p \ is the Planck mass and multiplying 1/M p i the 
scalar fields were rescaled to be dimensionless. The field- 
space metric for (j) 1 is assumed to be given by the 2x2 unit 
matrix. 

The ADM formalism has been utilized to derive the action 
of the dynamical variables particularly in the non-linear per- 
turbation theory ll43ll . Using the decomposed metric 

ds 2 = -N 2 dt 2 + hijidx 1 + 7VM<)(dx J + N j dt) , (2.2) 

the action is rewritten as 



$1* 



N S R- 2NV(ct> l , 2 ) + ^{E l3 E^ - E 2 ) 

(2.3) 



-Nh^difadjcp 1 



d A x 



where S R is the three-dimensional scalar curvature and Ei 
and E are defined by 

1 



{Pthij - DiNj - DjNt) , E h'-E, , (2.4) 



with the three-dimensional covariant derivative Di defined by 
hij. The spatial indices are raised and lowered by h%j . 

Perturbing the scalar fields as 0/ + <J0j, the background 
equations of motion are derived as 



fa + 3/50/ + dV/dtf = , 



6^ = ^ + 2^(01,02), 



(2.5) 
(2.6) 



where the dot denotes the derivative by the cosmological time. 
Using these equations, we also have 



9 = -^0/0 7 



(2.7) 



In this subsection, we describe our way to calculate the 
gauge-invariant quantity in the local universe. One simple 
way to preserve the gauge-invariance is to calculate fluctua- 
tions in a completely fixed slicing and threading. The time- 
slicing can be easily fixed by adapting the gauge condition 
at each spacetime point. By contrast, the complete fixing of 
the spatial coordinates is not easy-going in the local universe 
because of the degrees of freedoms in non-normalizable trans- 
formations [2H2O. We therefore seek for the way to construct 
the manifestly invariant quantity under the change in the spa- 
tial coordinates instead of fixing them completely. 

We fix the time slicing by eliminating the fluctuation in the 
adiabatic field. We also adapt the gauge conditions on the 
spatial coordinates, taking the spatial metric as 

hij = e 2 ^ + °[e s %, (2.8) 

where Sjij satisfies the transverse and traceless conditions 
5"f\ = diSYj = 0. While these conditions do not fix the 
spatial coordinates completely in the local universe, we need 
not worry about the presence of the residual gauge modes only 
if we consider the genuinely gauge-invariant quantity. We can 
construct the genuinely gauge-invariant operator by making 
use of the scalar quantities under the three dimensional spa- 
tial diffeomorphism. Following our previous works dHIH, we 
pick up the scalar curvature 3 R as such a scalar quantity. Al- 
though scalar quantities also vary under the gauge transforma- 
tion because of the change of their arguments x 1 , this gauge 
ambiguity does not appear in the n-point functions of these 
quantities with the arguments specified in a gauge-invariant 
manner. In order to specify the arguments of the n-point 
functions in a gauge invariant manner, we measure distances 
from an arbitrary reference point to the n vertices by means 
of the geodesic distance obtained by solving the spatial three- 
dimensional geodesic equation: 



dV 



dx 3 dx 
dA""dT 



0, 



(2.9) 



where is the Christoffel symbol with respect to the 

three dimensional spatial metric on a constant time hypersur- 
face and A is the affine parameter. We consider the three- 
dimensional geodesies whose affine parameter ranges from 
A = to 1 with the initial "velocity" given by 



d^(X,A) 



(1A 



-C(A=0) 



A=0 



-5 7 (A=0)/2 



jX j . (2.10) 



We identify a point in the geodesic normal coordinates X 1 
with the end point of the geodesic, x l (X , A = 1). Expanding 
the global coordinates x % by the geodesic normal coordinates 
X 1 as x % (X) =: X % + Sx l (X), a genuine gauge invariant 
variable can be given by 



9 R(X) 



R(t, x\X)) 
Sx 11 ■ ■ ■ 8x tn 



E 

n=0 



d x , n s R(t, X 1 



(2.11) 
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where we introduced the abbreviated notation X := {t, X}. 
Similar quantities to our genuine gauge-invariant variable 
were studied in Refs. I37l44lll . while the authors of Refs. ll37l - 
l40h seem to have the different opinions regarding what are 
actual observables. 

In the following section, we use the equality « which is 
valid only when we neglect the terms which do not yield IR 
divergences in the one-loop corrections to ( 9 R 9 R) . When we 
write down the Heisenberg operator in terms of the interac- 
tion picture fields for the adiabatic and entropy fields, this is 
equivalent to keep only terms without spatial/time derivatives 
and the terms which include only one interaction picture field 
with differentiation. In deriving the gauge-invariance condi- 
tion from the IR regularity of ( 9 R 9 R), we neglect the possibly 
divergent loop corrections whose loop integrals are composed 
only of the entropy field, because they are irrelevant to the 

IR 

gauge effects. Using the equality ps, the geodesic normal co- 
ordinates X % are related to the global coordinates x l as 



(X) 



IR 



S-y/2 



Using Eqs. ( fTTTT i, d214l) . and (l27T6t . we have 

9 i? 2 ™ d 2 ( 2 -i>d 2 x l d x ^, 
9 i? 3 « d 2 c 3 -(2d 2 (x*d x ^ 



(2.17) 
(2.18) 



Tpd 2 {X i d xi )Q 2 + -i) 2 d\X i d xi fi>. 



(2.19) 



III. MODELS WITH THE PURE ENTROPY FIELD 

In this and next sections, we investigate the predictions of 
the genuine gauge-invariant perturbation theory in two kinds 
of models. In general, the adiabatic and entropy fields can be 
coupled even in liner order. We first study a particular model 
where these two-fields are decoupled at linear order, deferring 
the study of coupled models to the next section. 



iX> 



(2.12) 



Abbreviating the unimportant prefactor, we simply denote 
the scalar curvature S R as 



S R 



IR 



- 2C [e-*T%C 



(2.13) 



The gravitational wave field 5%j can also generate the IR di- 
vergent with (djij Sjm ) ■ We can easily confirm that these con- 
tributions are generated only through the interaction between 
the adiabatic field and the gravitational wave field and that the 
entropy field does not produce another potentially divergent 
terms with the gravitational wave field. Therefore, repeat- 
ing the same argument as in Ref. [3|], we can show that the 
contributions from the gravitational waves are canceled in the 
gauge-invariant operator 9 R. To avoid the repetition of calcu- 
lations, we here do not explicitly write down the contributions 
from the gravitational wave field. 

Expanding C as C = Ci + C2 + C3 + ■•■ where we denote £1 
as tp := (1, the spatial curvature S R is expressed as 

s i?i = d 2 ip , S R 2 w d 2 ( 2 - 2i>d 2 ip , 

S R 3 « d 2 ( 3 - 2(C 2 d 2 V + ^d 2 ( 2 ) + 2i; 2 d 2 ^ . (2.14) 



A. Non-linear perturbations 

In this section, we consider the two-field model where the 
field 0i dominates the background. We use the horizon flow 
functions ll42ll defined in a similar way to those in single field 
models as 



2p 2 ' 



-m+l 



1 der m 



form>l. (3.1) 



Assuming that the horizon flow functions e m are all small of 
0(e), we neglect the terms of 0(e 3 ). We also employ the 
assumption: 



^ = 0(e 2 ). 
01 



(3.2) 



This condition requests that, within our approximation, the 
field <f>i and (j> 2 become the pure adiabatic and entropy fields, 
respectively. We define the gauge-invariant quantity 9 R on the 
slicing: 



5(j)x = 



(3.3) 



Using Eq. (12.12b . the difference between the geodesic normal 
coordinates and the global ones is expanded as 



5x l = Sx\ + Sx 2 + 



where 



5x[ w -^X l 



s.r: 



IR 



-C 2 x* + 



(2.15) 



(2.16) 



To evaluate the genuine gauge-invariant variable, it is con- 
venient to perform the calculation in the flat gauge: 



=2p IV7 



V 5 \ 



(3.4) 



where all the interaction vertexes are explicitly suppressed by 
the slow-roll parameters 01 EH. We associate a tilde with the 
metric perturbations in the flat gauge to discriminate those in 
the comoving gauge. Here again we can neglect the contribu- 
tions from the gravitational wave field, which do not appear in 
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9 R. The action in this gauge is given by 
r2 



N-\-6p 2 + 4/5<9 l iV l ) 



+ 



N- 1 ^ +ip T - N l d lVl ) (j> T + (p 1 - Wdup 1 ^ 

- 2N V yfl "/'> /l - ■ V m - Nh l: >d,wid 1 tp I 
ml 

(3.5) 

where we introduced V/ 1 .../ n := d n V/d(j> 11 ■ ■ ■ d<fr n and 
(50/ := ip/. Using the background equations and Eq. ( 13. 2k 
the derivatives of the potential are quantified as 



^2 



0(e 3 ), 



0( £ 2 ) . (3.6) 



The curvature perturbation in the gauge <50i = is related 
to the fluctuation of the dimensionless scalar field (divided by 
Mpi) in the flat gauge ipi as 

A IR /■ , /- a /■ , £ 2 . 2 , CACn 
C « Cn + Cn9 p C„ + —C n + j— 



where we have introduced 

Cn--=-(p/4>l)<pl- 0.8) 

Variation of the total action with respect to ipj yields 



~ 3p d t 



e 3p 



m=0 



7/7 



) - Ne- 2p d 2 ipi « . (3.9) 

Variations with respect to the lapse function and the shift vec- 
tor yield the Hamiltonian constraint: 



(n 2 -i)v+n 2 y 

' ml 

rn—l 

+ 2pd l N l + </>/^ 7 
and the momentum constraints: 



V 1 • • V" 



IR 



2^ 



(3.11) 



From the form of the constraint equations, we can confirm the 
presence of the degrees of freedom which appear from bound- 
ary conditions in solving these equations. Repeating a similar 
analysis to the one in Ref. J2l, these degrees of freedom are 
found to be the residual gauge degrees of freedom, which can 
change the average value of £„ at each order in perturbation. 
These constraint equations are solved to give 

5N g - ei6 , + | C 2 + \siS2 (C + *l) , (3.12) 



P 



1 



£l£2 (CndpCn + VndpUn) 



1 , , 2 

(2?722 + 3e 2 )cr„, 
where, in a similar manner to £„, we introduced 

cr„ := -{p/4>i) ip 2 ■ 



(3.13) 



(3.14) 



Substituting Eqs. (l3~T2l and (l3~T3l l into Eq. d3T9l , the evolu- 
tion equation of £„ is recast into a rather compact expression, 



L( n « -2e 1 C„ + -£i(4ei+e 2 )C 2 

-£l£ 2 Cn5 p Cn — ^£2£3Cn 7 

where the differential operator £ is defined by 

e -2p 



-2p 



£ := di + (3 - ei + £ 2 )d p 



(3.15) 



(3.16) 



Note that the adiabatic field is decoupled from the entropy 
field <j„ at linear order as the result of requiring Eq. (13. 2t . We 
also find that the interaction terms with £ n t7 n and £ 2<T " are 
suppressed at the order of 0(e 3 ). Since the entropy field s n is 
not included in Eqs. (13.7b and ( 13.15b , the curvature perturba- 
tion £ and the gauge-invariant spatial curvature 9 R are given 
in the same form as those in single field models. Therefore, 
in the decoupled model which satisfies Eq. (13.2b . the entropy 
field s n does not contribute to the IR divergence that origi- 
nates from the gauge effect. 

Reflecting these facts, the gauge-invariance conditions de- 
rived by requesting the regularity of loop corrections from the 
adiabatic field take the same form as those in the single field 
model. We briefly describe this result, deferring the detailed 
explanation to Ref. JH]. Here, we expand the interaction pic- 
ture field ijj as 



d 3 fc 



[ip k (X)a k +h.c] 



(2tt) 3 / 2 

with the creation and annihilation operators that satisfy 



(3.17) 



a k , a 



k' 



,5 (3) (fc-fc') 



(3.18) 



, (3.10) if me p 0S itive frequency function satisfy 



{1 + £l(l + £l + £ 2 )H - XiQ X* + £l + 2 £ 2 + 26 

- ^-L7 1 e 2 (2e 1 +e 2 )l MX) = -DMX) 



where Du and L k are defined as 
D k :=k-d k + 3/2 , 

L k := 2 + (3 - £i + e 2 )d p 



,-2p 



(3.19) 

(3.20) 
(3.21) 



the gauge invariant operator 9 R 2 is compactly summed as 



IR 



(l + X 2 )d 2 D k ^. 



(3.22) 
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Here £2 and A2 are time-dependent functions of 0(e 2 ) and 
they appeared as degrees of freedom in solving Eq. ( 13.15) at 
second order in perturbation. We can also find similar degrees 
of freedom at third order in perturbation and we label them by 
time dependent functions £3 and A3 of 0(e 2 ). If we choose 
them correctly fl, the gauge invariant operator 9 i?3 is com- 
pactly rewritten as 

IR 1 

m 3 « -iP 2 d 2 [(1 + 2\ 2 )D\1> - fiD k ^] - 5{ ni2 d 2 D k i> , 

(3.23) 

where we defined 5£ n ,2 and fx as 

Kn,2 :=-£- 1 ^e 2 (2ei+e 2 )V' 2 , (3.24) 

H ■= £1 + \ei - 3e? + -£i£ 2 + 2(£ 2 + A3) . (3.25) 

We finally find that the possibly divergent terms become the 
total derivative form as 

^R(X X \ W(X 2 )}) 

x{k 7 MXiWk(X2) + (cc.)} 

- (<5C«,a) / ^^d Xogk {k 7 MXiW k {X2) + (cc.)} , 

(3.26) 

and hence they vanish. Here we symmetrized about X\ and 
X 2 . We thus obtain the gauge-invariance conditions by re- 
questing the regularity of the loop corrections from the adia- 
batic field. Note that the preservation of the genuine gauge- 
invariance guarantees the regularity of the loop corrections 
from the adiabatic field, but it does not guarantee the regu- 
larity of those from the entropy field. Therefore, the divergent 
terms with (x 2 ) are still left in the left-hand side of Eq. (13.26) 
and they should be regularized by another method. 



B. Primordial non-Gaussianity 

In this subsection, we study the primordial non- 
Gaussianity, requesting the gauge-invariance conditions. In 
calculating the tree-level bi-spectrum, we can construct the 
genuine gauge-invariant variable using the curvature pertur- 
bation f in stead of the spatial scalar curvature S R. If we con- 
sider the modes with k 3> 1/ Lobs> where i b s denotes the 
observable scale in the comoving coordinates, the curvature 
perturbation £ evaluated in the geodesic normal coordinates: 

%(X) :=C(p, x\X)) 

= ((X) + Sx'diCU^x* +■■■ (3.27) 

becomes the gauge-invariant operator at least up to the second 
order in perturbation |]4|. 



Introducing the Fourier modes of 9 ( as 

we calculate the bispectrum of 9 £s at the leading order in per- 
turbation. Expanding % as 9 ( k = t/j k + 9 Q k 2 + ■ ■ • , we have 

(3.29) 

In the calculation of the bi-spectrum for 9 £, we again use the 

equality w, which picks up gauge-dependent terms. Using 
Eqs. ( 12.16b and ( 13.7b . the gauge-invariant curvature perturba- 
tion is expressed as 

% « Cn + Cn (d P - X'dxi ) C„ + \e 2 £ , (3.30) 

where ( n is given by solving Eq. ( 13.151 . Since there is no 
contribution from the entropy field in Eqs. ( 13.15) and ( 13. 30) . 
9 ( takes the same form as in single field models. As expected 
from these facts, repeating the same calculation as in Ref. fl, 
we again find that, in the squeezed limit with k\ <C k 2 ~ k^, 
the bi-spectrum vanishes as 

jp 

( 9 C fel WC fe3 > «0. (3.31) 

The bi-spectrum calculated in the conventional perturbation 
theory is thus found to be dominated by the gauge artifact 
also in the two-field model which satisfies the condition ( 13.21 ). 
Note that, while the three-point function generated from the 
cubic interaction £;] is suppressed from the request of the gen- 
uine gauge invariance, the one generated from the interaction 
(n&n i s not suppressed at all. Therefore, our gauge-invariance 
conditions do not affect on the three-point function with one 
adiabatic field and two entropy fields. 

One important remark is in order regarding the applicabil- 
ity of our argument. Precisely speaking, the expression of 5x l 
given in Eq. ( 12.16) can be reliably used only for the modes 
with k -C l/i bs- We can however show that the additional 
terms which appear for k > 1/ L Q b s vanish separately in the 
squeezed limit yfl. Thus, we can confirm that the result given 
in Eq. ( 13.31) is robust as long as we consider the super-horizon 
modes with k <C e p H. During inflation, the scales relevant to 
the current observations go far outside of the inflationary hori- 
zon, which leads to the natural assumption that the observable 
scale is much larger than the horizon scale 1/H. Therefore, 
our result is applicable to the observable modes which satisfy 
l/L obs < k < e p H. 

IV. GENERAL EXTENSION 

In the previous section, we studied the two-field model of 
inflation where the adiabatic and entropy fields decouple at 
linear order in perturbation. In this model, the allowable form 
of interactions is significantly restricted so that the entropy 
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field does not contribute to the potentially divergent terms 
with the loop integrals of the adiabatic field. The predictions 
of the genuinely gauge-invariant perturbation are then found 
to be essentially same as those in single field models of infla- 
tion. In this section, we consider more general two-field mod- 
els where the adiabatic and entropy fields are coupled even 
at linear order. The study of these models would help us to 
understand how the entropy field contributes to the IR correc- 
tions affected by non-local gauge degrees of freedom. 



A. Adiabatic and entropy perturbations 

In multi-field models, we need to care about the fact that the 
IR modes include the physical degrees of freedom as well as 
the gauge degrees of freedom. The discrimination of the later 
from the former proceeds easily in the model with the pure 
adiabatic and entropy fields, but in general it becomes more 
complicated. To ease this procedure, we decompose the fluc- 
tuations in the two fields into the horizontal and orthogonal 
directions to the background trajectory as 



5a 
6s 



-<t>2 4>1 



e 



8(p 2 



, (4.1) 



where we introduced the 2x2 time-dependent rotational ma- 
trix with a defined as 



(0l) 2 + (0 2 ) 2 



In this expression, the action (12. 3t is rewritten as 



N S R - 2NV(a, s) + — {!:'• I '., , - E 2 ) 



+ ^-{a + 5&- N t d l 5a - 65s) - Nh ij d i 5ad j 5a 
+ (6s - N 1 d l 6s + 68a) 2 - Nh ij diSsdjSs 



(4.2) 

where 6 is the local rotation angle given by 

8 —tan" 1 ^/^). 

The model studied in the previous section corresponds to the 
particular case with 8 = const, where the background trajec- 
tory is not curved. 

Using 6a and 5s, we rewrote the perturbation expansion of 
the potential as 

V(a, s) = V + V a 5a + V s Ss 

+ I^ + ^ s + i^ + ..., 

where Vq denotes the background value of the potential. Note 
that V ai .. ai .. an where = a, s are not the derivatives of 
the potential in terms of a and s, but they are given by the 
linear combinations of Vi 1 ,,i n . In general multi-field mod- 
els, the cross correlation between the adiabatic and entropy 



fields (SaSs) does not necessarily vanish. The loop correc- 
tions with (5a6s) can also yield divergent terms which are 
possibly affected by gauge degrees of freedom. Keeping this 
in mind, here we consider the terms that can contribute to the 
loop corrections with (5a 2 ) and also (5a5s). As in the previ- 
ous section, we keep terms that include at most one interaction 
picture field for 5a or 5s with differentiation, but among them 
we neglect terms that are expanded only in terms of the inter- 
action picture field for 5s. In the following we use the same 

IR 

equality rj to denote a equality which is valid in neglecting 
these terms. 

We define the gauge-invariant scalar curvature 9 R on the 
time-slice fixed by the gauge condition: 

5a = 0. 

We again change the gauge into the flat gauge given by 
Eq. d3.41 i. Noticing the fact that the background evolution 
is characterized only by the adiabatic field, we introduce the 
horizon-flow functions as 



2p 2 



-m+l 



1 dg m 
: m dp 



for m > 1. (4.3) 



In the following, we assume that the change in the background 
trajectory takes place satisfying d6 /dp = 0(e). We also as- 
sume that the mass of the entropy field satisfies V ss /p 2 = 
0(e) so that the entropy field participates in the generation of 
the fluctuation. In the flat slicing, the action is given by 

M 2 f 

S~—^ I dt d z xe* p N~\& + 5a - N l di5a - 96s) 2 

+ N-\5s - N'd.ds + 65a) 2 - Nh l] d l 5ad J 5a 

- Nh ij di6sdj6s + ^^(-Qp 2 + Apd^ 1 ) - 2NV . 

(4.4) 

The constraint equations are derived as 



2pd l iV i 



a5a + — (5a — 



1 



(5s + 66a) 2 



+ (N 2 - l)Vb + VJa 
1 



and 



-(V aG 5a 2 + 2V as 5a5s + V ss 6s 2 ) + 



2pdiN - N(a - 65s)di6a - N85adi5s 



IR 



0, (4.5) 



T TJ 

- 2N(5ad l 5a + Ssd.Ss) « , (4.6) 

where, in Eq. d4.51 l. we abbreviated the higher-order terms in 
V(a, s). In deriving the Hamiltonian constraint (14.5b . we used 
V s = —(7(9, which is given by varying the action with respect 
to 5s. Consulting Eqs. d4.5l ) and d4.61 l. we again find that the 
degrees of freedom in the boundary conditions are intrinsi- 
cally attributed to the gauge degrees of freedom in the adia- 
batic fluctuation, while the entropy fluctuations are, at non- 
linear order, also affected by the gauge modes through the 
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change in the lower-order adiabatic fluctuation. Varying the 
action with respect to 8a, we have the equation of motion for 

5a as 



z ~ 3p d t \ ^-(ct + 5a - 65s) \ - (a + 5a - (>5s) — 



+ N[V a + V tja 5a + V as 5s 

!t, r 2 1 

+ ^V aa<y 5a 2 + - 



Vaaa5a 2 + ^V ass 6s 2 + V aas 5a5s 



TD 

-~-(5i + 95a) - Ne~ 2p d 2 5a « . 
N 



(4.7) 



From a similar calculation to that in Appendix A of Ref . , 
the relation between the fluctuations in these two gauges is 
obtained as 

. ir , . „ j. 1 >9 d9 . 



1 



2 Wn + ^C^pCn + ^£2 (£2 + 2e 3 )£ 

3\d P ) ^n'W™ dp U " " 

-H^ + 2 £2 d7j c « s " + (d7) c " s "- (4 ' 8) 

where we defined Cn and s„ as 

Cn := -(fi/ar)Sa , s„ := -(p/&)Ss . (4.9) 

While in the case with d0/dp ^ the evolution equations 
for Cn and s„ are coupled even at linear order, we can easily 
solve Eq. (14.7b up to 0(e). In the following calculations, we 
keep the terms up to this order, neglecting the terms of 0(e 2 ). 
Then, the constraint equations ( 14. 5t and ( 14. 6t are solved to 
give 



~ TR 1 TR 

(SJV « -£lC„ , -ftJV* « £l9 p Cn 

P 



(4.10) 



Using Eqs. ( 14.1 Ob , the equation of motion for Cn is recast into 



ir AS. 
V 



d9 V a 



L( n « 2^<9 p S„ + ( 3^; - ^ ) S„ - 2e X Cn (?(n , 

(4.11) 



dp 



-2p 



where the derivative operator £ apparently takes the same 
form as the single field case: 



L := d 2 + (3 - El + e 2 )3 ( 



p -2 



(4.12) 



Up to this order, interaction terms with the entropy field do 
not appear in the equation of motion for Cn- Allowing the 
introduction of homogeneous solutions, Eq. ( 14. lit is solved 
to give 



IR 



+ f^ipX + M3"0 2 X + A2X 2 + AsV'X 2 ,(4.13) 



where £j, /x,, and /tj for i = 2, 3 are time-dependent functions 
of 0(e). Noticing the fact that the interaction picture field of 
s„ satisfies 

■Cx = 0(e), (4.14) 
we also included the homogeneous solutions with \. 

B. IR regularity and gauge-invariance conditions 

Now, we are ready to derive the gauge-invariance condition 
for the coupled case. Using Eqs. (12.18b and ( |2.19t together 
with Eqs. (gj) and J4.13t . we have 



m 2 « ,pd 2 



{(l + e 1 )d p -X i d xi 



26+£ 2 /2}v+(m2-^-)x 



— JV + 2/i 2 x 



M2 



dp 



(4.15) 



and 



^S 1 



-(d p -x*d xi ) 2 ip + 3&ip + » 3X 



(d p - X'dx^e^p + 36 + 3e 2 /4)^ 





d0\ 

M2 "d7> ; 



dfl 



M2 ~ ~d~p) ( " dp ~ X ' ldx ^ 
+ fi 2 (d p - AT l a x .)x + M3V' + fax 



(4.16) 



We keep the terms in the last line of Eq. ( 14.15b and the terms 
in the last two lines of Eq. ( 14.161 ), which can yield the possibly 
divergent terms with (ipx) ■ 

Making use of the Gram-Schmidt normalization, we can 
prepare a set of orthonormalized mode functions, which sat- 
isfy the Klein-Gordon normalization. We expand tp 1 := 
— (p/^Scf) 1 by the orthonormalized mode function ip 1 k 
where a = 1,2 are the indices for the orfhonormal bases. 
Using these orthonormal bases, tp and X are expanded as 



^(X) = Y J ®WJX): 
a=l 
2 



(4.17) 
(4.18) 



where the time-dependent matrix 9 is already given in 
Eq. d4.U and we defined ip^(X) as 



€(X) 



d 3 fc 

(2tt) 3 / 2 



(4.19) 
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with the creation and annihilation operators that satisfy 



®>a, k: Q>n k 



(4.20) 



Substituting Eqs. ( 14.17b and ( 14.181 ) into Eqs. ( 14.151 ) and ( 14.161 ). 

the gauge-invariant spatial curvature is expanded as 



m 2 « v Yl d ' 2 [ Q u{0- + ei)d p - X l d x , 

a=l 

+ 2&+e 2 /2} + e 2 jMal^ 



5Ia 2 [( M2 -— )ei/ + 2A 2 e 2/ ]^, (4.2i) 



and 



a=l 



-Q u {d p -X'd X i 



6i/(36 + 3e 2 /4 + ei9 p )(a p - X*fi^ ) 

e 2 /M3 p - x 4 9 x o + %9ij + M 3 e 2/ 

2 



r\f)\ 

e 1I [ f , 2 -—yd p -x i d x < 



a=l 

<d 2 l(l 2 (dp - X % d X i) + M3©1/ + A3©2/ 



(4.22) 

where we noted dOi//dp = (d9/dp)Q 2 i . 

In requesting the absence of the possibly divergent terms 
with (ifj 2 ) and (tpx), we obtain the gauge-invariance condi- 
tions as 



0ij{(1 + e x )d p - X l d x , + 2£ 2 + e 2 /2} + 821/12] < fc 
= -e u D k ii kl (4.23) 
£3 = M3 = , (4.24) 



and 



M2 



A0_ 
dp 



A2 = M3 = . 



(4.25) 



If all these conditions are satisfied, the gauge-invariant spatial 
curvature 9 R is simply given by 



w 2 ™-ijJ2d 2 D k e u ^ a , 



a = l 

2 



(4.26) 
(4.27) 



Then, the potentially divergent terms in 9 R are found to be- 
come the total derivative as 



({ fl fl(Xi), «R(X 2 )}> 

IR 1 

2 

xE{^ lfe (^i)<;ft) + (cc.)} , (4.28) 

a=l 

and they vanish. We now understand that, after the choice of 
the appropriate initial condition, the two-point function of the 
gauge-invariant operator 9 R is shown to be regular, leaving a 
possibly IR divergent contribution with (x 2 ), which is irrele- 
vant to the gauge effect. 

We requested Eq. (14.25b to eliminate the possibly divergent 
terms with (tpx)- However, at first glance it may be unclear 
whether the realization of the gauge invariance truly requests 
the absence of the terms with (tpx) or n °t, because the cross- 
correlation (tpx) i s a l so influenced by the entropy field. To 
make this point clear, we note that (tpx) no longer diverges if 
at least one of ip and x is suppressed in the IR limit. If we 
remove the residual gauge degrees of freedom, say by adapt- 
ing the local gauge condition [1], the adiabatic field is sup- 
pressed in the IR limit so that the regularity of {tpx) is suf- 
ficiently guaranteed. This indicates that the cross-correlation 
(ipx) does not diverge in the genuine gauge-invariant quanti- 
ties. Therefore, we can request the absence of the IR diver- 
gence from (ipx) as me necessary condition for the genuine 
gauge invariance. 

Several remarks are in order regarding the gauge-invariance 
conditions ( 14.23b - (14.25b . In harmony with the result in sin- 
gle field models, the gauge invariance condition (14.231 ) almost 
uniquely determine the mode function for ip to that for the 
Bunch-Davies vacuum at the leading order in the slow-roll 
approximation. This can be confirmed by following a simi- 
lar argument to the one in Refs. J2,|3t]. The gauge-invariance 
conditions thus derived should be consistent with the equation 
of motion for the adiabatic field. Using the Fourier expanded 
basis: 



1 



k 3 /' 



fi,k(p) e 



ih-X 



(4.29) 



the gauge-invariance condition ( 14.23b can be recast into 
{(l + ei)a p + fc-ft, + (2f 3 -ei)}e u /i = 0, (4.30) 



where we also used the first equation in Eq. (I4.25t . From 
Eq. ( 14.1 II ). we have the mode equation at liner order as 



(4.31) 



where we introduced the derivative operator L k as 

l k := ^ + 3(1 - ei)d p + e -^-k 2 - 3(£! + e 2 /2) 



(4.32) 
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Operating {(l+ei)d p +k-d k } onEq. ( 14.311 ). which commutes 
with L k up to 0(e), we obtain 

L k {(l + e 1 )d p + k-dk}QufLk 



where is defined as 



,d0 V as 



d6 



'% ~ ~f + 2 d~p dp ) {9p + k ' dk ^° 2lf ^ = ■ 

(4.33) 



The mode function for the entropy field \ is not restricted 
from the gauge-invariance condition. However, to reproduce 
the gauge-invariance condition (14.30b from Eq. (14.33) , we 
need to employ, at the leading order in the slow-roll approxi- 
mation, the condition 



(d p - X l d x , + D k )Q 2I ^ aM = 0(e) 



(4.34) 



which almost determines the mode function for x to that for 
the Bunch-Davies vacuum. Then, the terms in the second line 
of Eq. (14.33b vanish, reproducing Eq. ( 14.30b after the multi- 
plication of L^ 1 ■ The last two terms in ( 14.30b appear as the 
homogeneous solutions of L k . Following the same argument 
as in Ref. iH, we can also show that the gauge-invariance con- 
ditions and Eq. (14.34b sufficiently ensure that the mode equa- 
tions are satisfied for all wavenumbers, if they are satisfied 
only a particular wavenumber. 



C. Primordial non-Gaussianity 

Now we calculate the bi-spectrum for 9 £, which again be- 
comes the genuinely gauge-invariant variable for 1 /L ^ s <§; 
k <C e p H, in requesting the gauge-invariance conditions 
(I4T231 - (|4~25T >. Substituting Eqs. ( TTTB , d4~8t , and (gTTJi into 
Eq. (13.27b . the gauge-invariant curvature perturbation 9 £ is 
given by 

% « V + V>(1 + £i)d P 4> - ^x*d x ^ + Qe 2 + V , 

(4.35) 

where we again abbreviated the several terms in Sx l which can 
give the non- vanishing contributions for k ^> l/i b s . This is 
because we can show that these contributions independently 
vanish in the squeezed limit, repeating the same argument as 
in Ref. fl] . Note that after the use of Eq. ( 14.25b . the terms with 
the entropy field disappear in the expression of 9 (. The expres- 
sion of g ( includes the time-dependent function £2, which is 
restricted from the request of the consistent quantization JH. 
Assuming £2 is determined appropriately, we do not discuss 
the explicit form of £2, because this is not important for our 
discussions. 

Using Eqs. (14.17b and (14.18b . the gauge-invariant curvature 
perturbation 9 ( can be expanded as 



2 

a=l 



(l+e l )d p -X l d x , + (ie 2 + fc 



(4.37) 



The bispectrum for % is again given by the three terms 
in Eq. ( 13.291 ). where only the first two terms give dominant 
contributions in the squeezed limit k\ <C k 2 ~ £3. Using 
Eq. ( 14.361 ). the first term in Eq. ( 13.29) is given by 



IB 



QuQu «fe 1 ^,fe 2 M «'^ I fe 3 ) , (4-38) 

a,P,a',0'=l 



where M Q ^j,fc is the Fourier mode of JA a p- Note that 
(V'q fej ^ fc 2 ^a'/3',fc 3 ) gives the non-vanishing contributions 
only if a and j3 agree with either a' or (3', while a and /3 do 
not necessarily coincide with each other. We first consider 
the case with a = /?, where the non-vanishing contribution is 
given by 



in 



d 3 p 



u aM u a,k 2 



3 -i(fci+fc 3 )-X 



d 3 X 



(2tt) 3 / 2 J (2tt) 3 



X 



+ {2v)-*/H^\K)vi k y aM r aM d P r a , kl , (4.39) 

where ipa^ is defined as ip a . k '■= Qnv^ k and K is defined 
as K := ki + k 2 + k$. In deriving this expression, we used 
the gauge-invariance condition ( 14.23b and d p v^ = 0(e) for 
k <C e p H . Next, for the case with a ^ /?, we obtain 

(■^l.k^PM ( M "/3,fc 3 +M /3ct ,fc 3 )) 

S ,j v j f d3 p f d3x 



<*M U PM] ( 27r )3/2_/ ( 27r )3 

- i(hMX 6^\k 2 +p)r aM D P r^ x 

+ e-^ +k ^ x 5^\k 1+ p)r pM X l d x ^l p e lp - X 

+ (2Tty*iH^\K)vi k y^ M d P r aM . (4.40) 



a, 0=1 



The calculation of the second term in Eq. (13.29b proceeds 
similarly to give 

(V>fei s Cfe 2 ,2V>fe 3 ) 

2 

« QuQij J2 « k^p'M^ik 3 ) > (4.41) 
where the non-vanishing contributions are related to those in 
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the first term as 



V. CONCLUSION 



Hi 



61/61.7(7/4 fc 1 ^ ifc3 M QQ ,fc 2 ) 
- (2nr 3 ^S^(K) \i, aM \ 2 ^ aM d p r a , kl ]*(4*2) 

and for 

©l/©IJ<^a,ki (Mattel +^/3a,fc a )^,fc 3 ) 



IR 



(27r)- 3 / 2 <5( 3 )(K)|^, fe3 | 2 ^, fci a p ^ 



(4.43) 



Now, we consider the squeezed limit ki -C — ■ After 
some manipulation, Eq. ( 14.391 ) multiplied by Gi/0u can be 
expressed as 

^ IV'a.fcil 2 V'a,fc 2 C|fc 1 +fc 3 | 

x (27r)- 3 / 2 (fci-fc 2 + fc 3 )-^ 3 )(^) 
- |Va, fel | 2 |^, fc2 | 2 (27r)- 3 / 2 fci -d K 8^(K) 
+ (2tt) ~ 3/ V 3) (K) ^ a , fcl |Va,fe 2 | 2 P < fcl ■ (4.44) 

(The detailed calculation can be consulted in Ref. 0].) The 
terms with fci manifestly vanish in the squeezed limit k\ — > 0. 
The rest terms on the first and second lines are not sup- 
pressed at this moment. Combining the contribution from 
the second term of Eq. ( 13.291 l, which satisfies Eq. ( 14.42b . 
these terms, however, provide the factor (V>a.fc 2 '0a |fei+fe 3 | — 
"0Q,|fei+fc 2 | V'a fc 3 )> which again vanishes in the squeezed limit, 
fci — > 0. We therefore understand that only the term on the 
last line in Eq. ( 14.44b yields the non-vanishing contribution in 
this limit. Repeating a similar calculation, we again find that 
all the terms in Eq. ( |4.40t except for the last one vanish in 
this limit. Finally, the bispectrum for % in the squeezed limit 
fci —> is given by 



m 2(27r)" 3/2 ,5 (3) (K) 



2 

E 

a, 13=1 



|V'a,fe 2 | 2 Re[^ ( 3 jfci a p ^,fc 1 ] , 



(4.45) 



Note that for d0/dp = the right-hand side of Eq. (14.45b 
vanishes, reproducing the result in Eq. (13.31) , because for this 
particular case ip a .k becomes constant in time at superhorizon 
scales. By contrast, for d9/dp 7^ the adiabatic field does 
not necessarily become constant. The result in Eq. (14. 45b in- 
dicates that the time variation in the local rotation angle can 
generate the observable fluctuation, which are not eliminated 
by gauge transformations. 



In this paper we studied the implications of the genuine 
gauge invariance in two-field models of inflation. We showed 
that, likewise in single field models, if initial quantum states 
satisfy the gauge-invariance conditions, the loop diagrams 
with the adiabatic field do not yield IR divergences. It is 
remarkable that the gauge-invariance conditions, imposed on 
the adiabatic field, can be influenced by the entropy field. This 
is because the interactions between the adiabatic and entropy 
fields can generate the additional possibly divergent terms, 
which are absent in single field models, and these contribu- 
tions are also influenced by gauge effects. For the derivation 
of the gauge-invariance conditions, we distinguished the IR 
divergences which are relevant to gauge effects from those 
which are irrelevant to gauge effects, considering the dia- 
grams up to one-loop order. This discrimination would be- 
come rather complicated if we extend our argument to higher 
order in loops. The loop correction of the entropy field (8s 8s) 
then comes to be no longer gauge invariant due to the contam- 
ination of the adiabatic field. 

In requesting the gauge invariance in the local universe, we 
reexamined the bispectrum for the primordial curvature per- 
turbation. The conventionally used curvature perturbation £ 
preserves the invariance under normalizable gauge transfor- 
mations, but it does not under non-normalizable gauge trans- 
formations. This indicates that the curvature perturbation £ 
does not preserve the gauge-invariance in the local universe 
while this should be preserved in observable fluctuations. We 
therefore calculated the tree-level bi-spectrum for the gen- 
uine gauge-invariant curvature % to discuss observable fluc- 
tuations. In contrast to the result in single field models, where 
the genuine gauge-invariant bi-spectrum completely vanishes 
in the squeezed limit, in multi-field models, we still have the 
non-vanishing contributions in this limit. This is generated 
from the time variation in the curvature perturbation at super- 
horizon scales, which occurs only if the background trajectory 
is curved yielding d6 /dp 7^ 0. This effect can be understood 
as follows. While the constant part of the curvature perturba- 
tion can be removed by a local dilatation, which is the non- 
normalizable gauge transformation, the time dependent part 
cannot be removed by gauge transformations and produce the 
physical effect. We should emphasize that our arguments can 
be applied to the fluctuations in the observable scales of cur- 
rent measurements. Since the gauge-invariance conditions de- 
rived here are the necessity conditions, precisely speaking, the 
genuine gauge invariance of the bi-spectrum for % has not 
sufficiently proven. It is, however, intuitively reasonable to 
expect the non-vanishing physical effects in the presence of 
the entropy field. In this paper, we employed the slow-roll 
approximation. It would be interesting to discuss the observ- 
able effects in models which break the slow-roll approxima- 
tion, where large non-Gaussianities are predicted in the con- 
ventional perturbation theory. 

At the end of this paper, we add several comments on the IR 
divergence from the entropy field. The loop corrections of the 
entropy field can be divergent if the entropy field has the scale- 
invariant or red-tilted spectrum. This divergence cannot be 
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eliminated also in genuine gauge-invariant quantities. The IR 
divergence from the entropy field should be regularized by a 
different way from the one for the adiabatic field. In Ref. 113611 . 
we proposed one way to regularize the IR divergence from 
the entropy fields. (Some other ways of regularization were 
discussed in Refs. JUBMl ) We showed that, if we take into 
account the effects of the quantum decoherence which pick 
up a unique history of the universe from various possibilities 
contained in initial quantum state, the IR loop corrections of 
the entropy field no longer diverge. Therefore, if we consider 
the decoherence effect, it would be possible to show the IR 
regularity of the genuine gauge-invariant quantities. 
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